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We demonstrate the existence of new nonclassical corre- 
lations in the radiation of two atoms, which are coherently 
driven by a continuous laser source. The photon-photon- 
correlations of the fluorescence light show a spatial interferene 
pattern not present in a classical treatment. A feature of the 
new phenomenon is, that bunched and antibunched light is 
emitted in different spatial directions. The calculations are 
performed analytically. It is pointed out, that the correla- 
tions are induced by state reduction due to the measurement 
process when the detection of the photons does not distin- 
guish between the atoms. It is interesting to note, that the 
phenomena show up even without any interatomic interac- 
tion. 

PACS numbers: 42.50.Dv, 42.50.Ar, 03.65.Bz 



Resonance fluorescence from a single atom driven by a 
coherent field was the first example allowing to observe 
nonclassical effects such as antibunching and subpoisso- 
nian statistics These experiments first done with 

atoms in a beam have later also been performed using 
single trapped ions . Squeezing in resonance fluores- 
cence has also been investigated Q . Since the early work 
of MoUow I?), Carmichael and Walls |^ one has stud- 
ied at length the quantum statistical characteristics of 
radiation produced by a cooperative system of two and 
more atoms [p|jlC|] . Some evidence of cooperative effects 
have been demonstrated . There are also some recent 
proposals to investigate cooperative effects including the 
interactions between the atoms [lT2|,f3l . 
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FIG. 1. Two atoms in a coherent laser field show interfer- 
ence fringes 



In the present paper we study the two atom fluores- 
cence leading to new nonclassical effects. We consider 
a situation, where the two atoms are driven coherently 
by a continuous laser source. These nonclassical effects 
are observed in photon-photon-correlations when state 
reduction occurs in the measurement process. We note, 
that in an ealier study by Mandel correlation effects 
in a similar system have been discussed. However, in con- 
trast to our paper there no continuously pumped atoms 
have been considered. We calculate nonclassical two pho- 
ton correlations for all times including dissipation of the 
atoms. The dynamics resulting from the excitation by 
cw fields and by spontaneous emission is very important 
for the present results. We should mention that on the 
experimental side first order interference in the radiation 
produced by a prepared system of two atoms has been 
reported |l^. We also point out, that the creation of 
entangled states of distant atoms has been discussed in 
a different context JlGf . 

It may be noted, that the trapped ion technology is 
fairly well advanced to envisage resonance fluorescence 
measurements with chains of trapped ions driven by co- 
herent fields. One can also use single molecules in crys- 
talline hosts, where recently nonclassical photon statis- 
tics have been measured ]l7[ |. 

We consider two identical atoms with the levels \e)i and 
\g)i {i = 1, 2) at fixed positions xi and X2 with dipole mo- 
ment d and transition frequency m. They are driven by a 
resonant external laser field with wavevector k/,. We as- 
sume that the only dissipative terms are due to the spon- 
taneous decays of the levels |e)i,2. In the rotating-wave, 
Markov and Born approximations, the time evolution of 
the system is given by the master equation |18 



p^-\^^Y^ (kL ■ x„ - <^ii)] ^+ ^ 



c.,p 



~ X] ^ (^m'^M P + /"^M " 2(7^ pcr+) , (1) 

where 2r2 is the Rabi frequency of the atom laser sys- 
tem, (7^ are the atomic raising and lowering operators 

for atom number /i, and 27 = ||dp^ is the Einstein A 
coefficient for the single atom. As we assume the elec- 
tromagnetic field outside the laser beam in the vacuum 
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state, the correlation functions can be written as 
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GW(r,i)= e[''=^--^l«(t)a;(t)) (2) 

G(2)(ri,t;r2,i + r)= ^ Jifc(x,,. ■ ri + xp.^ ■ r^)] ^ 

(a+(i)<(t + r)a;(i + r)a;(t)) (3) 

with x^^^ := - Xi. and f = j7| ■ 
The two time correlations (H) can be calculated using 
the quantum regression theorem and the time evolution 
of the density operator. Therefore the normalized second 
order correlation function can be written as 



J2) 



(ri,t;r2,t + r) = 



P(r2,^ + T|ri,t) 
P{r2,t) 



(4) 



where P(r, t) is the probability of finding a photon at 
position r at time t, and P{r2,t + T\ri,t) the conditional 
probability of finding a photon at position r2 at time t + r 
premising the detection of a photon at point ri at time 
t. Using the inequality 

Tr(p(aI(ri)+/3I(r2))')>0 Va,/3 (5) 



leads us to 
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(6) 



for classical systems. This gives us a good possibility 
to estimate the nonclassical behavior of our system. 

As there can always be found one state which does not 
interact with the laser field it is suggestive to use this 
state to build up the basis for the further calculations. 
We call this state la) and define 



|e) := |e,e) 
V2 



IS) := — [e ■•^|e,g; + e"^|g,e; 



\a) :=-L(e-</'|e,g)_ei<^|g,e; 
|5> := |5,5> 



with ■■=\i)i®\i)2, 0:=-kL-xi2 



(7) 




FIG. 2. Level scheme of the atomic system for the sym- 
metrized set of states 



In this representation the master equation reduces to 
the following sets of 6+3 equations 

Pee = 4(apg^ - Pee) 

Pss = 2(pee - Pss + 2a(p!,g - p^J) 
Paa = 2(pee - Paa) 

PIs = -3PL - 2a(pee - Pss + Pig) 

pig = 2p^^ - p\g + 2a(l - Pee - Paa " 2pss + Peg) 

Plg^~2{plg + a{plg- p"^,)) (8) 

Pla = -3Pea " 2apL 

p'sa = 2(a(p::, + p^g) - p: J 

Pag = -2p^„ - 2ap^„ - pig 



(9) 



with a:^^, pli := 5R(pfe;) and p^^^; := 3(pfe;)- 



By solving the first set of equations (||) we find the 
diagonal elements in the steady-state to be: 
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ee 


(72 + 2^2)2 



(10) 

(11) 

PeT = . . . _..2 (12) 



This calculation does not depend on the direction of 
the driving laser, although it is included by the proper 
definition of our symmetric and antisymmetric state. 




FIG. 3. Steady state populations of the symmetric and an- 
tisymmetric state as a function of the laser intensity Q 



If the laser direction is perpendicular to the ion separa- 
tion (0 = 0) this solution corresponds to the earlier found 
solutions of Richter where we neglect the dipole- 
dipole interaction. This is justified when the ion dis- 
tance is in the range of several wavelengths of the atomic 
transition. In this case the solutions are completely inde- 
pendent of the separation of the atoms. For strong laser 
fields (O 7) the populations are equal for each of the 
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atomic states [p^^ 



ss 



P, 



ss 



ss 

raa 



1/4). 



We next discuss, how a detection event leads to state 
reduction and entanglement. After a detection at the 
point r we have to find a new density operator in the 
following way 



p(r) 
a±(r) 



(^+(r)T-(r)) ' 



where 



(13) 
(14) 



If, for instance, the system is in state |e) before a de- 
tection at r with fxi2 = 0, we find it in \s) after the 
detection. This explains, how entanglement could come 
into the system, as the system is transformed from a 
non-entangled state |e) to an entangled state \s) by a de- 
tection event. 



In this case the system behaves as the well known dou- 
ble slit experiment where the atoms act as slits and the 
two optical paths interfere. Like in a system with pre- 
pared initial states the photons emitted by symmetric 
transitions |e) — > |s) and |s) \g) show inverse fringes 
compared to those emitted by antisymmetric transitions 
|e) — > \a) or \a) \g). This follows from the 1 ± cos(5(r) 
terms in (p^. One major difference is that the contrast 
depends on the intensity of the laser light so that the 
fringes disappear for higher laser intensities when f2 ^ 7. 
The populations pff and pf^ equalize. 

To get the second order correlation function we solve 
the master equation by calculating the Liouville Opera- 
tor and use for the initial state the density matrix given 
by (|l^). This procedure is also equivalent to using the 
quantum regression theorem. 



Let us now check, if the atoms remain uncorrelated 
after a detection event i.e. p{r) factorizes into a product 
of density matrices pi with pi referring to the density 
matrix of the i-th atom. We find that this is not the 
case, for example 



^{p{r)sa) = 



n'^ sin ( 



2(f^2 + ^2(i + cos(5(r))) 



,Pee(r)-0 (15) 



with S{r) = (ki — fc^f) • X12. On the other hand a 
factorized density matrix would imply, S(p(r)so) = af- 
ter any detection event. Furthermore, the dependence of 
Q{p{r)sa) on the coordinates of the two atoms through 
(5(r) also points out the nonexistence of a factorized den- 
sity matrix after detection. Thus there are values of r, 
where the state reduction by the measurement leaves the 
system in an entangled state. With this mechanism the 
entanglement, which is necessary for nonclassical effects, 
is brought into the system simply by the detection of a 
single photon. 

To calculate the intensity at position r we use the re- 
lations 

{afa^) = = p^e + ^ {Pss + Paa) (16) 

(a+a^-) = {a+a^Y = \^'t> (p., - p,,) , (17) 

where we assumed the system to be in the steady state 
(i.e. pf^ = 0). We then find 

G(i)(r,t) = (l + cos<5(r)) [p,S) + Pss{t)] + 

{I - cos 5{y)) [p,,{t) + paa{t)]. (18) 

This leads to the well known interference fringes in the 
first order correlation function as studied in the experi- 
ment of Eichmann et al. [|5| . 



Remarkably enough we are able to give an analytical 
expression for the intensity-intensity correlation: 
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5(2)(ri,0;r2,r) = 1 + w , y\ ^ 

Av^{s + cosdi)(s + COSO2) 

{4e^*i/^ssinJi sin^2 + s [e*i^^s + (s - 1)^] cosJi cos (52 

-e^^/'^vs^ [2iycos(i/t) +3sin(zyt)] 

-|-2e^*/^:/s (cos(5i -|-cos(52) [(2s- 3) sin(z^t) - 2j/cos(i/t)] 
-f e*/^ V (2e* cos 5i cos 82 + s sin 5i sin 82) 
[2v{s - 2) cos{vt) + (5s - 6) s\n{vt)] 



+ - cos 61 cos 82 [(s(s(4s - 33) + 64) - 36) cos{2vt) 
+2v{s - 2)(5s - 6) sva{2vt)] ] 



(19) 



with s := 2 +1^5^ = 5(r,) and v iV8s- 9 

being parameters for the intensity of the driving field and 
t := 7r is the time scaled with the atoms spontaneous 
emission rate. For r = Eq. ( |is| ) reduces to 



7(2)(ri,0;r2,0)- 



s2 cos2(^ 



(s -I- cos 5i){s + cos(52) 



(20) 



Assuming measurements with only one detector at po- 
sition r, ^'^2) (j-^ 0; r, 0) shows interferences fringes with 
maxima larger and minima smaller than one (Fig Q). 
This means that the system is emitting super-poissonian 
light in one and sub-poissionian light in the other direc- 
tion. This behavior has no analogy in the first order 
correlation i.e. in the two slit experiment. 
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FIG. 4. Second order correlation function for one detector 
with Si = 82 = S, Q = 0.87 

For two detectors the situation is different. For any de- 
tector position with 62 = (2n + l)7r + Si, g^^-'(ri, 0; r2, 0) 
vanishes completely, while at other positions we find 
maximum detection probability (Fig ^). After a detec- 
tion at position ri with Si = 2m: we find pee = Paa = 0. 
Therefore there is no probability to detect any photon at 
a position with S2 = (2n -I- l)7r at the same time, as no 
emission into the antisymmetric channel could take place 
(p^). So we again find fringes as a function of ^2. For 
di = {2n + 1)tt we detect an emission on the antisymmet- 
ric channel first, so we find g^^)(ri, 0; r2, 0) vanishing for 
positions r2 with S2 — 2mi for the same reasons as above 
(Fig§. 




FIG. 5. Second order correlation function for two detectors 
with &i = 2n7r, = O.87 




FIG. 6. Second order correlation function for two detectors 
with 5i = (2n + l)7r, Q. = O.87 

After any detection on a specific channel the system 
needs time for re-excitation to emit to the orthogonal 
channel; on the other hand cascade emissions could only 
take place at the same channel. Note that the behav- 
ior of (7*^^)(ri, 0; Y2, r) is well understood in terms of our 
symmetrized basis (0). These results show also that the 
inequality (^) could be violated by choosing the two 
detector positions in a way that g^'^\ri, 0; ri, 0) > 1 and 
f;'^^^(r2, 0; r2, 0) < 1. In this manner we find nonclassical 
behavior in our system as consequence of the detection 
induced entanglement of the two atoms. We note that 
dependence of the type cos^ in the two photon correla- 
tions has also been derived by Mandel in his work on 
radiation from atoms prepared in atomic coherent states 
[Q. In contrast to that we are dealing here with the 
more general case of continuous excitation of the atoms 
including dissipation. 

In conclusion it has been shown by analytic calcula- 
tions, that a system of two atoms coherently excited by 
a cw laser source shows nonclassical features being de- 
pendent on the distance of the atoms, although there 
is no assumed interaction between the atoms. This can 
be understood by measurement induced entanglement, 
which comes into the system if the detection does not 
distinguish between the two atoms, i.e. the detected pho- 
ton does not carry which-way-information. On the first 
glance the system studied seems to be similar to the two 
slit experiment, however the nonclassical effects discussed 
here do not show up in the first order correlation and can 
only be seen in the second order correlation as it should 
be the case for a quantum phenomenon. 

To see the reported effects single ions stored in a lin- 
ear rf-trap can be used; remaining micromotion could 
be overcome by a phase-sensitive detection In this 
paper the case was discussed, that the detector/or de- 
tectors do not discriminate the light emitted by the in- 
dividual atoms. There are, however, further interesting 
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phenomena observable when the which-way-information 
is available. One example is that the excitation of one 
atom and the selective observation of the fluorescence 
from the other one opens, among other effects, the possi- 
bility to investigate the dipole-dipole interaction between 
the atoms. These results will be described elsewhere. 
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